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Abstract 

We discuss the mass modifications of D and D* (B and B*) mesons in nuclear medium. The 
heavy quark symmetry for D and D* (B and B*) mesons is adopted, and the interaction between a 
D or D* (B or B*) meson and a nucleon is supplied from the pion exchange. We find the negative 
mass shifts for D meson and B meson, and hence that the D and B mesons are bound in the nuclear 
medium. As applications, we consider the atomic nuclei with a D meson, ^°Ca and g 08 Pb, and 
investigate the energy levels of the D meson in each nucleus. We also discuss the mass shifts in 
the isospin asymmetric nuclear medium, and present the possible phenomenon about distribution 
of isospin density around a D or B meson in nuclear medium. We find that the mass shifts of D* 
and B* mesons have large imaginary parts, which would prevent precise study of the energy levels 
of D* and B* mesons in nuclei. 

PACS numbers: 12.39.Fe,12.39.Hg,14.40.Lb,14.40.Nd,21.65.Jk 



I. INTRODUCTION 



Research of hadrons with charm and bottom flavors is one of the most interesting subjects 
in the present hadron and nuclear physics. As motivated by discoveries of new exotic hadrons 
with charm and bottom (X, Y, Z etc.) in experiments, new types of hadrons, such as multi- 
quark or hadronic molecules, are investigated by many researchers p~H5]. Among them, 
hadrons including D or B meson are interesting objects not only in vacuum but also in 
nuclear medium. Because a D (B) meson is composed with anti-charm quark c (anti-bottom 
quark b) and a light quark q = u, d, there is no annihilation process of quark and antiquark 
pair in nuclear medium. The decay modes are only electromagnetic or weak processes. 
Therefore, we expect to obtain precise information about the dynamics of D (B) meson in 
nuclear matter, such as their energy levels in atomic nuclei, without being disturbed by 
the strong decay processes and absorption processes [6H8]. This nice property of D and B 
mesons may be in contrast with their antiparticle states, D and B mesons. Because D and 
B mesons in nuclear medium include light quark annihilation processes, the nonnegligible 
widths by the decays and absorptions make it difficult to study the dynamics of D and B 
mesons in nuclear medium. 

Study of D and B mesons in nuclear matter will provide us with important information 
about (i) hadron-nucleon interaction, (ii) modification of hadrons in nuclear matter, and 
(iii) changes of nuclear medium caused by hadron as impurity. As known in the light flavor 
mesons in nuclear medium, the mass shifts of vector mesons u, p, and give information of 
partial restoration of chiral symmetry breaking in nuclear matter [9] . It is also the case for D 
and B mesons, because each of them contains a single light quark. The properties of D and 
B mesons are also concerned with the modifications of quarkonia in nuclear medium [TOl - IT^] . 
In early works [SHE], it was discussed that D and B mesons are bound in atomic nuclei, such 
as 208 Pb. The properties of D mesons have been studied by several researchers with various 
theoretical approaches, such as QCD sum rule jTUHPT] , hadronic dynamics [T5H23] and so on. 
In the literature, the binding energy of a D meson in normal nuclear matter is estimated to 
be around a few ten MeV. In the present paper, we discuss the hadron dynamics by focusing 
on two important symmetries for heavy-light mesons; the heavy quark symmetry for a heavy 
quark and chiral symmetry for a light quark. 

For the heavy quark symmetry, as a general property of QCD, the magnetic gluon couples 
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to the heavy quark with a suppression factor 1 / itiq with the heavy quark mass tjiq . It means 
that the spin of the heavy quark changes with a suppression by order of 1/wt.q, and hence 
it leads to the approximate mass degeneracy of pseudoscalar meson and vector meson in 
hadron mass spectrum [26H30] . The mass splitting between D and D* is around 140 MeV. 
The mass splitting between B and B meson is around 45 MeV, hence the heavy quark 
symmetry becomes better for bottom flavor. Indeed, these mass splittings are smaller than 
those in light flavor mesons; about 400 MeV between K and K* mesons, and 630 MeV 
between 7r and p mesons. Therefore, due to the mass degeneracy in heavy-light mesons, 
we need to take int account both pseudoscalar mesons and vector mesons simultaneously as 
effective degrees of freedom. The importance of the approximate degeneracy of D and D* 
(B and B*) mesons was already emphasized in studies of several exotic charm and bottom 
hadrons, such as two-body DN (BN) systems [3TW34"] . Zb mesons [35], T cc states [36] [59] 
and so on. It will be shown in the text that the approximate mass degeneracy is important 
also for the dynamics of D and D* (B and B*) mesons in nuclear medium. 

Another important symmetry is chiral symmetry. The pion dynamics induced from the 
spontaneous chiral symmetry breaking in vacuum provides us with a key to understand 
the low energy hadron and nuclear physics |37j . In recent studies, it is discussed that the 
nuclear systems are almost governed by the pion dynamics following to the strong tensor 
force. For example, the binding energies of nuclei are given about 80 % by the pion exchange 
as shown in the few-body calculations [3H]- The nuclear matter with infinite volume can 
also be described almost by the pion exchange between nucleons [39-44]. Not only the 
normal nuclei, but also neutron-rich nuclei [45] and hypernuclei with strangeness j4"6l44"8"] . 
are important objects in which the pion dynamics plays a dominant role. The importance of 
the pion exchange was also emphasized for exotic charm and bottom hadrons in Refs. [31f - 
[36] as above mentioned. Especially, the various bound or resonant states of the two-body 
DN (BN) systems can be formed by the strong tensor force induced from the pion exchange 
between D (B) and N [31TI34] . Therefore, it is quite essential to consider the pion dynamics 
in order to discuss the properties of D and D* (B and B*) mesons in nuclear medium. 

In the present work, with respecting the heavy quark symmetry and chiral symmetry, we 
focus on the pion dynamics for D and D* (B and B*) mesons in nuclear medium. It will 
provide us with information about (i) D- and D*-nucleon (B- and B*-nucleon) interaction 
and (ii) modification of D and D* (B and B*) mesons in nuclear matter. We discuss shortly 
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also (iii) changes of nuclear medium caused by D and D* (B and B*) mesons as impurity. 

This paper is organized as follows. In Sec. 2, by following the description of the pion dy- 
namics in nuclear matter [3"5H4"Tj and hypernuclear matter [IBHIS] , we apply their approach 
to D and D* (B and B*) mesons in nuclear matter. We adopt the heavy quark symmetry for 
DD*7r and D*D*7r (BB*7r and B*B*7r) vertices, and obtain the self-energies of those particles 
in isospin symmetric nuclear matter. We find that D and B mesons have negative mass 
shifts at the normal nuclear matter density, and hence they can be bound in nuclear matter. 
In Sec. 3, we apply the results in the uniform nuclear matter to the study of the energy 
levels of D meson in atomic nuclei. We consider g°Ca and g 08 Pb nuclei, in which a D meson 
is embedded in Ca and 208 Pb nuclei. We also discuss isospin asymmetric nuclear matter 
and find the upper and lower components of the isospin doublet of D (B) meson behave 
differently. We discuss the possible phenomena about the distribution of isospin density 
around D (B) meson in nuclear medium. In Sec. 4, we summarize our discussion and give 
some perspectives. 

II. D AND D* MESONS IN NUCLEAR MEDIUM 

D and D* (B and B*) mesons containing a heavy quark with a light antiquark have two 
different important symmetries; the heavy quark symmetry for a heavy quark and chiral 
symmetry for a light quark. Both symmetries are incorporated in the effective theory for 
heavy-light mesons [2"oT - [3T)] . The interaction Lagrangian with pion is given as 

= -i TrH a v,iV£ a H b + g n TrH a H bW Al, (1) 

where we define the heavy-light meson field H a via pseudoscalar field P a (D or B meson) 
and vector field P* M (D* or B* meson), 

with the four dimensional vector v with constraint v 2 = 1, isospin index a, and the Dirac 
indices \l — 0, • ■ • ,3. The vector and axial vector currents are V 1 = |(<^<9 M £ + £c^£ ) and 
= |(£t#*f - £#*ft) W ith f = e iM/f and 



The coupling constant g n = 0.59 is determined from the DD*7r decay width. / = 132 MeV 
is the pion decay constant. The above vertex form is rewritten as 

A™* = j^Pv u (7r-d u 7r + i(n x d^)-f)P^ - j^P*^(ir-d u n + i(n x d v T?)-r)P$ (4) 

+y^p;t + h.c. - 2 ^p : ^p;P d »M ba e aXPv v\ 

For nucleon parts, we define the NN7r vertex as 

^ = -Jf2^(« x d ^>? N + fjN W d^-rN, (5) 

with the coupling constant #a = 1-3. Then, the interaction between D or D* and nucleon is 
given by the pion exchange [60]. The present formalism will be applied to B and B* with 
better accuracy thanks to their heavy masses. 

Let us consider the self-energies of D and D* mesons in uniform nuclear matter at zero 
temperature. The leading contribution is given by the two pion exchange of which relevant 
diagrams are presented in Figs, [I] and [2] [61J. However, the diagrams in Fig. [2] give zero 
contribution because of the cancellation of momenta at vertices or the antisymmetrization 
of isospin indices in isospin symmetric nuclear matter. Therefore, we consider only the dia- 
grams in Fig. [TJ We may additionally consider the terms in which the two pion propagators 
are crossed. However, such diagrams are suppressed by the inverse of the nucleon mass or 
heavy meson mass. We neglect such higher contribution in the present discussion. In the 
nucleon loop, we use the in-medium propagator for nucleon with mass ton and momentum 
p = (po,p) in isospin symmetric nuclear matter with Fermi momentum fcp, 

W + m N ) I-, V— - " 27r <V " 00(Po)0(k F -\p\)\ l f , (6) 

in which i] is an infinitely small positive number and If is a 2 x 2 unit matrix in isospin 
space for proton and neutron [3"9TT4T| HHH1H]- The second term in Eq. ^ is represented 
by the horizontal double line in Fig. [TJ We assume, for simplicity, that the uniformity of 
nuclear matter does not change by the existence of the D and D* mesons. The discussion 
beyond this simplicity will be given later. The nucleon mass in the propagator is set to be 
equal to the mass in vacuum. The imaginary part is not adopted in the in-medium nucleon 
propagator at zero temperature. The Fermi momentum hp is related to the nuclear matter 
density n = 2kp/3n 2 in isospin symmetric nuclear matter. 
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We give the expression of the self-energies for D and D* mesons. The same forms are 
applied to B and B* mesons by replacing their masses [62J . The self-energy of D meson is 



E D (A*) = £ D (M^ C *) + E D (A*)<£*T\ (7) 

with 



\(27rB*) _ Q&Sa f d ^ f d Vi rnsW 



/ 4 J (27T) 3 ./^l^ (27r)3 A 2 + U 1 + rm)\£\2 + ^.£_ irl (|£|2 + m 2)2 



from la) in Fig. [T] and 

/ -%2|<*F ( 27F ) ^|ft[<*F ( 2?r ) 



tmnIpi — P2I 



1 



from lb). The self-energy of D* meson is 
with 

S D ,(fc F )^ C *) = / ^ / ^ , ^ , (11) 

" ' (2^) 3 7 bl |< fcF (2VT) 3 l(l + m)|£|2 + ^. £ _^ ( | £ |2 +m 2)2 

d 3 p*i r« N |£| 4 1 

Pl|< fcF (2vr) 3 _ A 2 + 1(1 + ^)\£\2 iv (|£| 2 + m2)2 ' 




(12) 



from 2a) in the same figure, and 



y. (h _^9a [ d % f d 3 ^ 

/ ./|fc[<*F ( 27F ) i|pi|<fc F ( 27r ) 

^N|Pl-P2| 4 1 



|(1 + =f)|p 2 | 2 + I(-l + ^Ipll 2 - ^p r p 2 _ foj - p 2 | 2 + m 2 )2 



^D^fcpjp u - / — — - / — — - (14) 

/ ^|p 2 |<fcp ( 27F ) i|pi|<fcp ( 27r ) 



^nIpi — P2I 



i 



-A 2 + |(1 + ^)|p 2 | 2 + 1(_1 + «M)\fr\2 _ ngfr.fr _ ir] Qfr _fr\2 + m 2)2 

from 2b). We note that there is no D state in the intermediate state in the self-energy 
of D meson, because the DE)7r vertex is forbidden due to parity conservation. We define 



A 2 /mN = mg. -mg as a mass difference between D and D* mesons, and M = {m,Q+3mQ*)/4 
as an averaged meson mass. In the integrals, £, pi, and f>2 are three-dimensional momenta 



first accomplish the four- dimensional integral and then expand by the inverse of the nucleon 
mass and D and D* meson masses. We consider only the leading terms in the expansion. 
It is known that the approximation works not only in normal nuclear matter but also in 
hypernuclear matter (39HH1 H6TH8] . It will work well also for charm and bottom systems 
because the masses are heavier than those of nucleons or hyperons. 

Because the integrals by £ diverges, the regularization should be performed. According 
to the procedure in Refs. (IH1 HTj, we adopt the following regularization for the real parts; 



with x = 0, ±A 2 . Here, P stands for the principal integration. The sum of the first and 
second integrals becomes finite in the integration at \£\ — > oo, while the last term becomes 
finite by introducing three-dimensional momentum cutoff parameter A. Thus, we succeed 
to separate the cutoff-independent (first and second) term and -dependent (third) term as 
discussed in Refs. [3"9"H4"Tl |4"61[4"8"] . The three-dimensional momentum cutoff parameters are 
set to be Ag = 1.27Aa and Ab = 1.22Aa for D and B mesons, respectively, with Aa = 700 
MeV for A hyperon in nuclear matter jl6]. The ratios of Aj/Aa or Ab/Aa are estimated 
by assuming that the cutoff momentum is proportional to the inverse of the size of hadron; 
Ap/Aa = ta/Yp with mean-matter-radii rp (r\) of P = D, B meson (A hyperon) which 
is determined from the quark model [31J. This relation is satisfied well for the ratio of A 
and nucleon. Indeed, we find that A n /Aa = 650 MeV/700 MeV = 0.93 to reproduce the 
properties of (hyper) nuclear matter [39TPTI1 06HIE] is comparable with ta/tn = 0.94 from 
the mean-matter-radius estimated in the quark model. In the numerical calculations, we 
use the following mass values; = 938.9 MeV, = 138.0 MeV, mg = 1867.2 MeV, 
m * = 2008.6 MeV, m B = 5279.3 MeV, and m B * = 5325.1 MeV. 

Before discussing the numerical results, we remark that the imaginary part of the D (B) 



for pion (£) and nucleon {p 1 , p 2 = p\ + £) in the loops. To obtain the above equations, we 




(15) 
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meson self-energy is zero at low density. Because there is a mass gap between D and D* 
(B and B*) mesons in the intermediate states in la,b) in Fig. [I] The imaginary parts are 
absent for n < n CT = 2(kp) 2 /3n 2 , where the critical Fermi momentum is given by 



(16) 



This can be checked, because the functions in the momentum integrals for the D (B) meson 
self-energy have no pole for k-p < k F v . The imaginary part appears only at high density, 
n > n cv . As numerical values, we obtain the critical densities n cr = 2.2 fm -3 for D meson 
and n CI = 0.28 fm -3 for B meson. In both cases, the critical densities are larger than the 
normal nuclear matter density n = 0.17 fm -3 . Because the present formalism cannot be 
applied at high densities, we do not need to consider such high density state and ignore the 
imaginary part of the D (B) meson self-energy. On the other hand, the imaginary parts in 
the self-energies of D* and B* mesons exist even at low densities. The value of the imaginary 
part is independent of the cutoff parameter. 



(1a) D 
D* 
D 

(2a) d* 
D*,D 
D* 



n 



n 



n 



(1b) D 
D' 
D 4 



n 



(2b) □* 

nD D *' D 
D* 



ID" 



FIG. 1: The diagrams of the self-energies for D (la, b) and D* (2a, b) mesons. The horizontal 
double lines represent the second term in Eq. Q. See the text for the details. 



We present the numerical results for D and D* (B and B*) mesons in the isospin symmetric 
nuclear matter. The self-energies of D and D* (B and B*) mesons are plotted as functions 
of nuclear number density in Figs. [3] and |4j respectively. We obtain —35.1 MeV for D meson 
and —106.9 MeV for B meson at the normal nuclear matter density n®. As mentioned before, 
there is no imaginary part. On the other hand, we obtain —153.9 — zT70.1 MeV for D* meson 
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FIG. 2: Another diagrams of the self-energy for (1) D meson and (2) D* meson. 

and —203.3 — H17.5 MeV for B* meson at the normal nuclear matter density n . Therefore, 
these states are unstable in nuclear medium due to the large decay widths, though they are 
long-lived in vacuum with small decay widths (less than a few MeV for D* meson and zero 
for B* meson in strong interaction). 

The negative mass shifts of D and D* (B and B*) mesons indicate that they can be bound 
in nuclear matter. This result will be reasonable, because the interaction between a D (B) 
meson and a nucleon N can be attractive by the pion exchange as discussed in Refs. [3"1TI31E] . 
There, not only the two-body DN (BN) systems, but also the possible three-body DNN 
(BNN) systems have been discussed (31]. Therefore, we expect that the D (B) meson and 
nucleon systems have a rich variety from small baryon numbers to infinitely large baryon 
numbers. 

III. DISCUSSIONS 

A. Energy levels of a D meson in atomic nuclei 

We consider the atomic nuclei containing a D meson. The potential of D meson in the 
nucleus can be given as the real part of the self-energy [63]. We assume that the properties 
(density distribution, shape, and so on) of the core nucleus do not change by the existence 
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FIG. 3: The real part of the self-energies for D and B mesons. 
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FIG. 4: The real and imaginary parts of the self-energies for D* and B* mesons, 
of D meson and adopt the Woods-Saxson-type potential for D meson 

v (r,A) = l + J'_ RA)/ „ (17) 

with nucleus radius Ra = r^A 1 ^ with baryon number A. The parameters are tq = 1.27 fm, 
a = 0.67 fm from the empirical data for stable nuclei (see for example Ref. |5T|). Vq is the 
potential at the center of the nucleus. From the results of the self-energies in the previous 
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section, we use Vq = —35.1 MeV for D meson. We solve numerically the Klein-Gordon 
equation in the non-relativistic approximation and obtain the energy levels of D meson in 
g°Ca and g 08 Pb, in which D meson is embedded in 40 Ca and 208 Pb nuclei, respectively [64]. 
As shown in Fig. [5J several s- and p-wave orbitals of D meson inside the nucleus exist. The 
binding energies for the ground states (g.s.) of D meson are 28.2 MeV for g°Ca and 32.8 
MeV for g 08 Pb. The latter binding energy is close to the absolute value of the self-energy 
-35.1 MeV of D meson in nuclear matter with infinite volume. 

40Q a 208p b 



2.0 MeV 1.9 MeV 

5.9 MeV 



14.9 MeV 



p-wave 9i2 MeV 



18.7 MeV 



7.3 MeV 
p-wave 



27.0 MeV 

28.2 MeV g 

s-wave 

32.8 MeV 
s-wave 



FIG. 5: The energy levels of D meson in ^°Ca and ^ 08 Pb. 



The energy levels for D* meson in atomic nuclei may be discussed similarly. In this case, 
we need to introduce the complex potential with a finite imaginary part. The complex 
potential can be regarded as the self-energy of D* meson in the previous section. However, 
the imaginary part is so large (~ 100 MeV) that the fine energy splittings (~ a few MeV) 
between the orbitals of D* meson in nuclei will not be seen clearly. In comparison with 
D mesons, therefore, D* mesons are not so useful objects for the study of energy levels in 
nuclei. 

So far we have considered D mesons. Similar discussion may be applied to their antipar- 
ticle states, namely D mesons. In fact the only modification is just a change of the sign of 
the vertices of DD*7r and D*D*7r. However, there appear additional strong decay processes 
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(DN — > 7rE c , vrS*, and so on) as well as two-body absorption processes (DNN — > A C N, S C N, 
£*N, A*N, and so on). Those channels must make the analysis be complicated and give 
finite imaginary parts of the self-energy of D meson [52]. The similar situation will be also 
the case for B meson. Therefore, we claim again that D and B mesons are unique hadrons 
to study the fine energy splittings of D and B meson orbitals in nuclei. 



B. Isospin asymmetric nuclear matter 

In isospin asymmetric nuclear matter, the Fermi surface is unbalanced due to the different 
number densities of protons and neutrons. We introduce the Fermi momenta k F and k F 
for protons and neutrons, respectively. Then, the step function in the in-medium nucleon 
propagator in Eq. (|6| is modified as 

6(k F - \p\)l f -> 6(k F - \p |)i±* + 6(k F - \p (18) 

where If is a 2 x 2 unit matrix in the isospin space, and (1 ± t%)/2 are projection operators 
for protons and neutrons, respectively [3D]. It is convenient to introduce the asymmetry 
parameter 5 by giving k F and k F as 

kr = k F (l T 6) 1/3 , (19) 

where 5 is equal to (n n — n p )/(n n + n p ) with proton or neutron number density n p>n . 5 = 
is the case of isospin symmetric nuclear matter, and 5 = 1 is the case of neutron matter. 
The total nuclear matter density is n = n n + n p with n = 2/c|/37r 2 and n P)n = (/cp' n ) 3 /37r 2 . 

When the nuclear matter is isospin asymmetric (5^0), we observe mass splitting of up 
and down components in isospin doublet of D = (D°, D~) or B = (B + , B°). We remark that 
the mass splitting of the isospin doublet is zero in the isospin symmetric nuclear matter. 
We obtain the self-energies of the upper and lower components in the isospin doublet D 



or B by using the in-medium nucleon propagator in Eq. (18) [65]. In Fig. 6 we present 



the self-energies of D° and D~ mesons as functions of the asymmetry parameter 5 at the 
normal nuclear matter density n p + n n = n = 0.17 fm~ 3 . We find that the mass splitting 
in the isospin doublet becomes larger as the asymmetry parameter increases. Interestingly, 
the self-energies of the upper and lower components in the isospin doublet of D meson are 
negative at 5 = 1, and they are still bound in neutron matter. We may notice that the mass 
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splitting increases linearly as a function of the asymmetry parameter 8 when \8\ is small. 
We can parametrize the self-energies of D° and as 

£g ym W) * Z { i ym \n) - a»(n)5, (20) 
E^\n,8) ~ E^ m \n) + a Q (n)8, (21) 

for small \8\ <C 1, where £g ym ^(n) is the self-energy of D in isospin symmetric nuclear matter 
and ±0p (n) 5 denotes the deviation. We obtain crg(n) as a function of nuclear matter density 
n in Fig. [TJ We find erg (no) = 30 MeV at the normal nuclear matter density uq. 

In Fig. pi we plot the results for B + and B° mesons also. For the B meson case, the 



Fermi energy of neutrons or protons reaches at the critical Fermi momentum in Eq. (16) at 
\S\ = 0.63, and the self-energies of B + and B° mesons have imaginary parts for \S\ > 0.63. 
For small \8\, we define a^{n) for B meson by giving the self-energies of B + and B° mesons 
as 

^7 m) (n,5) ~ Sj ym) (n) - a B {n)8, (22) 
Ej ym) (n, 8) ~ Sg ym) (n) + a B (n)8, (23) 

where Sg ym ^ (n) is the self-energy of B meson in isospin symmetric nuclear matter. We show 
(Jb(ti) in Fig. [7j and find ctb(jio) = 97 MeV at the normal nuclear matter density tiq. 

We notice again that the modifications of the self-energies of D and B mesons are pro- 
portional to 8. This property can be contrasted with the change of the binding energy of a 
nucleon; the change of the binding energy per a nucleon is proportional to 8 2 [30]. Hence, 
for small deviation of isospin asymmetry (\8\ <C 1), the energy gain of the upper (lower) 
component of D or B meson for 8 > (8 < 0) overcomes the energy loss of the binding 
energy of a nucleon. 

The above result leads to an interesting phenomena about the spatial distribution of 
isospin density around D or B meson in nuclear matter. Let us suppose that D or B meson 
with isospin I z = +1/2 (the upper component; D° or B + ) is embedded in nuclear matter. 
Here, we assume that the nuclear matter is uniformly isospin symmetric (8 = 0) before D 
or B meson is embedded. From a view point of energy minimization, it will occur that the 
isospin density becomes neutron-rich (8 > 0) in the region near D° or B + meson. Because, 
for small \8\, the energy gain by D° or B + meson dominates against the energy loss of nuclear 
matter. In order to conserve the total isospin in the whole nuclear matter, the region far 

13 



apart from D° or B + meson will become proton-rich, because that the region will be less 
affected by the existence of D° or B + meson. A similar phenomena would occur for D or B 
meson with isospin J z = —1/2 (the lower component; D~ or B°), when the role of isospin 
up and down is just interchanged. The spatial distribution of isospin density around D or 
B meson embedded in uniform nuclear matter may be called "isospin polarization." Thus, 
the isospin of D or B meson in nuclear matter would be dressed by the deformation of the 
isospin density by protons and neutrons. We present the schematic picture in Fig. [8j The 
self-energies in the previous section were obtained under the condition of the uniformity of 
the isospin density in the nuclear matter around D or B meson. We expect that, if the 
variation of the isospin density in nuclear matter is taken into account, the binding energy 
of D or B meson would become larger. However, a quantitative discussion will be left for 
future studies. 

The picture presented here may provide us with a new dynamics of D or B meson in 
nuclear medium. One can expect in general that the phenomena of "isospin polarization" 
may occur if the hadron embedded in nuclear medium as impurity has the following two 
properties; (i) absence of strong decay and (ii) isospin non-singlet. As for the first condition, 
we should note that D and B mesons would have large decay widths which prevent to form 
stable states like isospin polarization in nuclear medium. As for the second condition, we 
may consider in the strangeness sector that K or K meson (isodoublet) can be embedded 
in nuclear medium [53J. However, it is known that systems of K meson and nucleon are 
unstable because of the strong decays with large width and the nonnegligible absorption 
processes, although K meson may be deeply bound in nuclear medium [66]. K meson is 
not bound in nuclear medium. E, S c , and £b baryons (isotriplet) will be also unstable in 
nuclear medium, because they have strong decays by n emission and absorption processes 
with nucleons. Therefore, we conclude that D and B mesons are unique impurities which 
can exhibit the stable isospin polarization in nuclear medium. 

IV. SUMMARY 

We discuss the mass modifications of D and D* (B and B*) mesons in nuclear medium. 
With respecting the heavy quark symmetry and chiral symmetry, we calculate the self- 
energies of D and D* (B and B*) mesons by the pion exchange. As results, we obtain the 
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FIG. 6: The mass splittings of D = (D°,D _ ) (top) and B = (B+,B°) (bottom) at the normal 
nuclear matter density no = 0.17 fm~ 3 as functions of the asymmetry parameter 5. 

self-energies —35.1 MeV for D meson and —106.9 MeV for B meson in isospin symmetric 
nuclear medium at the normal nuclear matter density 0.17 fm~ 3 . Thus, both D and B 
mesons can be bound in nuclear matter. As an application to atomic nuclei, we find several 
energy levels of D meson in g°Ca and g 08 Pb nuclei, in which a D meson is bound in 40 Ca 
and 208 Pb nuclei. We obtain the self-energies —153.9 — H70.1 MeV for D* meson and 
—203.3 — H17.5 MeV for B* meson at the normal nuclear matter density. However, due to 
their large imaginary parts, the individual energy levels of D* and B* in atomic nuclei will 
not be seen clearly. We further discuss the isospin asymmetric nuclear matter, and find that 
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FIG. 8: Schematic picture of isospin density distribution around D° (I z = +1/2; left) and D 
(J z = -1/2; right). 

the upper and lower components in the isospin doublet D or B meson split in mass for any 
small isospin asymmetry. With this result, we present the possible phenomena of "isospin 
polarization" induced by D or B meson in nuclear matter. 

For future studies, we will need to investigate the following several effects, which are 
omitted in the present formalism. The discussions about the higher order effects of the pion 
exchange, the sigma term, A excitations, and so on, will be required for an extension to 
higher densities jUJ [33]. It will be also necessary to consider the modification of the QCD 
vacuum at finite density. In experiments the nuclei with D or B meson may be produced in 
collisions of hadron (antiproton, n) to nucleus |55] , in nucleus-nucleus collisions [56J and 
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in relativistic heavy ion collisions [571 58J. The future experimental studies in accelerator 
facilities, such as J-PARC, FAIR at GSI, RHIC, LHC, and so on, will be expected to be 
performed. 
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